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A theory for the equilibrium low-temperature magnetization M of a diluted Heisenberg antifer- 
romagnetic chain is presented. Only the nearest-neighbor exchange interaction is included, and 
the distribution of the magnetic ions is assumed to be random. Values of the magnetic fields B t 
at the magnetization steps (MST's) from finite chains with 2 to 5 spins (pairs, triplets, quartets, 
and quintets) are given for chains composed of spins S—5/2. The magnitudes of these MST's as 
a function of the fraction, x, of cations that are magnetic are given for any S. An expression for 
the apparent saturation value of M is derived. The magnetization curve, M versus B, is calculated 
using the exact contributions of finite chains with 1 to 5 spins, and the "rise and ramp approxima- 
tion" for longer chains. An expression for the low-temperature saturation magnetic field B s (n) of 
a finite chain with n spins is given. Some non-equilibrium effects that occur in a rapidly changing 
B, are also considered. Some of these result from the absence of thermal equilibrium within the 
sample itself, whereas others are caused by the absence of thermal equilibrium between the sample 
and its environment (e.g., liquid-helium bath). Specific non-equilibrium models based on earlier 
treatments of the phonon bottleneck, and of spin flips associated with cross relaxation and with 
level crossings (anticrossings), are discussed. Magnetization data on powders of TMMC diluted 
with cadmium [i.e., (CH 3 ) 4 NMn x Cdi_ x Cl 3 , with 0.16 < x < 0.50] were measured at 0.55 K in 18 T 
superconducting magnets. The field B\ at the first MST from pairs is used to determine the NN 
exchange constant J. This J/kB changes from —5.9 K to —6.5 K as x increases from 0.16 to 0.50. 
The magnetization curves obtained in the superconducting magnets are compared with simulations 
based on the equilibrium theory. A reasonably good agreement is found. Data for the differential 
susceptibility, dM/dB, were taken in pulsed magnetic fields (7.4 ms duration) up to 50 T. The 
powder samples were in direct contact with a 1.5 K liquid-helium bath. Non-equilibrium effects, 
which became more severe as x decreased, were observed. For x=0.50 the non-equilibrium effects are 
tentatively interpreted using the "Inadequate Heat Flow Scenario," developed earlier in connection 
with the phonon bottleneck problem. The more severe non-equilibrium effects for x=0.16 and 0.22 
are tentatively attributed to cross-relaxation, and to crossings (more accurately, anticrossings) of 
energy levels, including those of excited states. For x=0.16 (lowest x), no MST's were observed 
above 20 T, which is attributed to a very slow spin relaxation for pairs, compared to a millisecond. 
A definitive interpretation of this and some other non-equilibrium effects is still lacking. 

PACS numbers: 75.50.Ee, 71.70.Gm, 75.10.Jm, 75.60.Ej 



I. INTRODUCTION 



Spin clusters with predominantly antiferromagnetic 
(AF) interactions exhibit steps in the equilibrium mag- 
netization as a function of magnetic field. These magne- 
tization steps (MST's) arise from energy- level crossings 
which change the ground state. They are observed at 
very low temperatures when only the ground state con- 
tributes to the magnetization M . In recent years MST's 
have yielded a wealth of information about AF clusters, 
first in diluted magnetic materials and later in molecular 
magnetism. An overall review of MST's was published 
recently^ For recent reviews of the magnetic properties 
of molecular clusters, including MST's, see Ref. y. 



In a molecular crystal the AF clusters are normally all 
of one type. The MST's then give values of exchange con- 
stants and anisotropy parameters for that cluster type. 
A diluted magnetic material, on the other hand, contains 
numerous types of spin clusters. Different cluster types 
give rise to different series of MST's. In addition to ex- 
change constants and anisotropy parameters, the MST's 
also give information concerning the populations of the 
different cluster types. The populations are related to 
the magnitudes of the MST's in the different series. The 
results for the cluster populations can be used to check 
if the distribution of the magnetic ions is random. 

Most previous studies of MST's in diluted magnetic 
crystals were on three-dimensional (3D) materials*! al- 
though some quantum wells were also studied^ The 
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present paper, however, is devoted to MST's from a di- 
luted AF Heisenberg chain (ID). The material studied is 
TMMC [chemical formula: (CH 3 ) 4 NMnCl 3 ] which was 
diluted by replacing a large fraction of the Mn atoms by 
Cd. Powder samples of (CH 3 ) 4 NMn x Cdi_ x Cl 3 , with x 
between 0.16 and 0.5, were investigated. In these ma- 
terials the non-magnetic Cd 2+ ions break the chains of 
Mn 2+ ions into finite segments. 

Pure TMMC is probably the closest approximation to 
an ideal isotropic (Heisenberg) linear AF chain. For re- 
views of its magnetic propert ies, w ith extensive references 
to original works, see Refs.l4l l5l6l l7l The ID magnetic be- 
havior of this compound is due to the crystallographic 
structure. It contains chains of Mn ions which at room 
temperature are along the c-axis of the hexagonal struc- 
ture. The Mn ions in each chain are linked by CI ions. 
The space between the combined Mn — Cl 3 — Mn chains 
is occupied by tetramethylammonium groups. 

The strongest magnetic interaction in TMMC is the 
isotropic exchange between nearest neighbor (NN) Mn 2+ 
ions in the chain. The NN exchange constant, ob- 
tained from various experiments (e.g., Hef. Si. is J/fce = 
—6.6 K, where &b is the Boltzmann constant. Other 
intrachain exchange constants are believed to be much 
smaller, and are usually neglected. The exchange inter- 
action between different chains is orders of magnitude 
smaller than the intrachain interaction, but is responsible 
for the long-range AF order below the Neel temperature 
T w =0.84 K. The anisotropy in TMMC is mainly due 
to the dipole-dipole interaction. It is of the easy-plane 
type, and is two orders of magnitude smaller than the 
dominant exchange interaction^ 

At room temperature both TMMC and its Cd ana- 
log [(CH 3 ) 4 NCdCl 3 , known as TMCC] have isomorphous 
hexagonal structures (space group P6 3 / m ). Crystallo- 
graphic phase transitions at lower temperatures result in 
a lower symmetry, and in small structural differences be- 
tween TMMC and TMCCi?' 1 ^ 11 ' 12 These differences are 
often assumed to be unimportant, although it is conceiv- 
able that even small changes in the crystal structure have 
some effect on the magnetic behavior, especially on spin 
relaxation at low temperatures. 

Previous investigations of (CH 3 )4NMn x Cdi_ x Cl 3 in- 
cluded measurements of the susceptibility and the Neel 
temperature as a function of x (Ref. tlSh . These re- 
sults were interpreted theoretically both by the origi- 
nal authors, Dupas and Renard, and by Harada el alM> 
Susceptibility measurements on a related diluted linear 
AF chain (DMMC:Cd) were obtained and interpreted by 
Schouten el alM The authors of both Refs. [Tj and [Tj 
noted the difficulties of preparing alloys with uniform Cd 
concentrations . 



II. EQUILIBRIUM THEORY 

In this section a theory for the low temperature mag- 
netization M of a diluted AF Heisenberg chain is pre- 



sented. It is assumed that the spin system is in thermal 
equilibrium with a constant-temperature heat reservoir. 
This equilibrium theory is suitable for interpreting the 
data that were obtained in the slowly-varying magnetic 
fields ("dc fields") of the superconducting magnets. The 
additional considerations needed to interpret the results 
obtained in pulsed fields, of several ms duration, will be 
discussed in Sec. IIHI 



A. The Model 

MST's from finite AF chains were predicted decades 
agojiSiii but some more recent theoretical results are 
also useful. The simplest model for MST's in a diluted 
magnetic material is the single-./ cluster modeli It in- 
cludes only the largest isotropic exchange constant J and 
the Zeeman energy. Other exchange constants, and all 
anisotropics, are ignored. This model, with the NN in- 
trachain exchange constant chosen as J, is expected to be 
a good starting point for (CH 3 ) 4 NMn x Cdi_ x Cl 3 (here- 
after, TMMC:Cd). All cluster models are applicable only 
when x is not too high. However, for a diluted magnetic 
chain (ID) the single-J model it is expected to hold at 
least up to x=0.5. All the samples in the present study 
are in this range. 

In the single-J model the magnetic clusters are finite 
chains, each consisting of n coupled spins. These clusters 
are treated as independent. The total magnetization M 
is then the sum of the magnetizations of finite chains with 
different n. Let fi n be the average magnetic moment 
of a chain with n spins, and let N n be the number of 
"realizations," per kg, of a chain with n spins. (N n may 
also be called the "population," per kg, of finite chains 
with n spins). The magnetization per kg is then 

If P n is the probability that a spin is in a chain with 
n spins, and if iVtotal is the total number of spins per kg, 
then 

N n = N total P n /n. (2) 

The probabilities P n are obtained from well known 
results^ 

P„ = nx"- 1 (l-x) 2 . (3) 

This result assumes a random distribution of the mag- 
netic ions. This crucial assumption is discussed later. 
Figure H shows the probabilities P n for n < 5, and the 
probability P>5 that a spin is in a finite chain n > 5, i.e., 

5 

P>5 = 1 ~ = (6 - 5x)x 5 . (4) 

l 
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FIG. 1: The probability P n that a magnetic ion is in a finite 
linear chain with n spins. Results for n < 5 are plotted as a 
function of the fraction x of cations that are magnetic. Also 
plotted is the probability P>5 that a magnetic ion is in a 
cluster with more than 5 spins. 

From Eqs. (^3J), the magnetization M is given by the 
infinite sum 

M = AT total ^x n - 1 (l-x)V- (5) 

n 

In practice, exact results for the average magnetic mo- 
ment fj, n are available only for values of n which are not 
too large (short chains). The infinite sum in Eq. (JjJJ 
is therefore truncated at some maximum value of n, 
called 7i max . Chains with n > n max are treated using 
an approximation^ In the present work we chose n max =5 
because exact results for chains with up to five Mn 2+ ions 
(5=5/2) were readily available from previous worksii^ 

When x is below 0.50, less than 11% of the spins are 
in chains with n > 5. For this range of x, the "rise and 
ramp" approximationi can be used for the total contri- 
bution from chains with n > 5 [i.e., it approximates the 
remainder after the sum in Eq. J5| is truncated] . The ap- 
proximation smoothes the independence of this remain- 
der, i.e., MST's from chains with n > 5 are not resolved. 
This approach (exact treatment for n < 5, and an ap- 
proximation for n > 5) proved to be adequate for inter- 
preting the data obtained in dc magnetic fields. However, 
the interpretations of some data obtained in pulsed mag- 
netic fields also used available theoretical results^ for 
n > 5. 

To implement the rise and ramp approximation, the 
so-called "short chain model" (Ref. 15) was adopted. 
This model is appropriate for x < 0.5 because more than 
99.4% of the spins are in chains with n < 10, and more 
than 99.9% are in chains with n < 14. A chain with 14 
spins is still short enough to be described by that model. 
In the short-chain model the ground state at B=0 has 
total spin 5t(0)=0 when n is even, and 5y(0)=5 if n 



is odd. This simple result does not hold for 3D materi- 
als, which is one of the reasons why the rise and ramp 
approximation is simpler and much more accurate for di- 
luted chains than for diluted 3D materials. 

B. Qualitative shape of the magnetization curve at 
low temperatures 

1. Clusters of one type 

Consider first an ensemble of identical finite chains, 
all with the same n. (In the language of Ref. [lj, this is 
an ensemble of "realizations" of a chain "type" with n 
spins.) At low temperatures, fc^T <C |J|, the qualitative 
variation of \i n with B depends on whether n is even or 
odd. When n is odd the zero-field-ground-state total spin 
5t(0)=5 aligns rapidly at low B. This alignment is given 
by the Brillouin function (BF) for spin 5. The rapid rise 
of /i„ at low B ends when the BF approaches saturation. 
After the BF saturates, there is a magnetic-field interval 
in which fj, n is nearly constant. At still higher B a series 
of MST's appears. The fields at these MST's depend on 
n. Once this series of MST's is completed, fi n reaches its 
true saturation value ^„ imax =ng/iB5. 

When n is even, St (0)=0. Therefore, no initial fast 
rise of fi n occurs at low B, in contrast to the case of odd 
n. However, at high magnetic fields a series of MST's 
still appears. At the completion of this series, fi n reaches 
its true saturation value. 

2. Total magnetization M 

Chains with all values of n contribute to the total mag- 
netization M. When ksT «C | J|, the chains with odd val- 
ues of n produce a fast rise of M at low B. This rise fol- 
lows the BF for spin 5. After this fast rise is completed, 
and before the appearance of the first MST of significant 
size, there exists a field interval in which M stays ap- 
proximately constant. This (nearly) constant value, M s , 
is the "apparent saturation value" (see Ref. Q). At still 
higher fields, MST's series from chains with n > 2 ap- 
pear. Once the last MST of significant size is completed, 
M reaches its true saturation value Mq. As discussed 
later, the magnetic field required to saturate the magne- 
tization of a chain remains finite even when n — > oo, so 
that M reaches true saturation at a finite B. 

C. Apparent saturation value 

The true saturation value of the magnetization is 

M = AW15MB5. (6) 
The apparent saturation value is 

M s = ]T iVtotaix"- 1 (l-x) 2 . g/ iB5, (7) 

n— odd 



4 



where the sum is only over odd n. Therefore, 

M s /M = x"- x (l-x) 2 . (8) 



3(1(1 



This infinite geometric series can be summed^ 
M./M = (l-x)/(l + x). 



(9) 



It is noteworthy that an exact analytical expression 
for M s was not obtained for 3D materials. 1 A rough ap- 
proximation was then used for the net contribution from 
clusters with n > n max . This approximation was used 
only for x < 0.1. In contrast, Eq. 10 for a diluted chain 
is exact, within the framework of the short-chain model. 
It should be very accurate for x < 0.5. 



D. MST's from Chains with n < 5 

All chains with n > 2 give rise to MST's. The magnetic 
fields Bi at the MST's from chains composed of spins 
5=5/2 were given earlier for n=2 (pairs, or dimers), rt=3 
(triplets, or trimers), n=4 (quartets, or tetramers), and 
n=5 (quintets, or pentamers)^ For completeness, the 
values of the reduced fields 6, = g^Bi/\J\ are repeated 
here. 

For n=2 there are five MST's at 6^=2,4,6,8, 10. For 
n=3 there are five MST's are at 6 i= 7, 9, 11, 13, 15. For 
n=4 there are ten MST's at ^=0.95, 2.04, 3.39, 5.02, 
6.87, 8.85, 10.88, 12.94, 15.00, 17.07. For n=5 there are 
ten MST's at 6*=4.62, 5.89, 7.18, 8.49, 9.86, 11.29, 12.83, 
14.48, 16.24, 18.09. 

Figure [21 shows the zero-temperature values of fj, n as 
a function of b for n=l, 2, 3, 4, 5. At finite T the ground 
state of a chain is not the only contributor to /i„. It is 
then necessary to include all energy levels in the calcu- 
lations of fi n . The procedures for such calculations were 
discussed in Ref. 



E. Rise and ramp approximation for chains with 

n > 5 

The total contribution of chains with n > 5 to M is 
approximated by a sum of two terms: 1) a fast rise at 
low B, which follows the BF for spin 5, and 2) a linear 
"ramp" from B=0 up to an effective saturation field B s . 

The fast rise at low B is due to chains with odd n, 
starting with n—7. Its magnitude is 



(AM) risc = M x 6 (l-x)/(l|x). 



(10) 



The linear ramp approximates the superposition of nu- 
merous MST's from all chains with n > 5. The ramp is 
given by 



ill, 



ramp 



(AM) 
(AM) ra]:lj , 



P (B/B S ), for B < B s 
for B > B s , 



(11) 
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FIG. 2: Magnetic moment per cluster, fi n , for chains with 
n < 5 at T=0, plotted as a function of the reduced field 
b=gfi B B/\J\. 



where 



(AM), 



P >5 M - 
M x 5 (6 



(AAf) riS o 
-4x 2 )/(l+x). 



(12) 



The reduced field b s (n) where the magnetization of 
a finite chain with n spins reaches saturation at T=0 
increases with n. However, in the limit n — > oo it is still 
finitejiiiS namely, b s (n—oo) = 85. It can be shown that 
the n-dependence of b s {n) is given by the equation 



b a (n) = 8Scos 2 (iT/2n). 



(13) 



For chains with n > 5 the change of this function is only 
7%. In the present case of 5=5/2, b s (n) changes from 
18.7 to 20.0 when n increases from 6 to oo. The value 
6 S =19, corresponding to 



gti B B a = 19|J|, 



(14) 



will be used in Eqs. (|llfl for the ramp. 

Figure^a) shows the predicted zero-temperature mag- 
netization curve for 5=5/2 when x=0.50. Figure|^b) is 
an expanded view for the range of magnetic fields rel- 
evant to the present experiments. The integers in this 
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FIG. 3: a) Calculated zero-temperature equilibrium magne- 
tization of a diluted linear chain with x=0.5 as a function of 
the reduced field b. These results are for 5=5/2. b) Expanded 
view of part (a) for the range of reduced fields relevant to the 
present work. The values of n for the chains responsible for 
each MST are indicated. 



III. NON-EQUILIBRIUM EFFECTS 
A. Non-Equilibrium Effects 

Experiments are normally conducted with the sample 
in contact with, or near, a thermal reservoir of constant 
temperature, e.g., a liquid helium bath at a tempera- 
ture Tbath- In some experiments the sample is not in 
thermal equilibrium either within itself and/or with the 
thermal reservoir. That is, the time for reaching com- 
plete equilibrium (both internal and with the reservoir) is 
not short compared to the time of the experiment. Such 
non-equilibrium cases require special considerations. 

Thermal equilibrium is usually maintained if 1) the 
magnetic field B is swept slowly ("dc fields"), and 2) 
the sample is in good thermal contact with the reservoir. 
These conditions are often fulfilled in superconducting 
magnets when the sample is immersed in liquid helium. 
There are, however, exceptional cases of non-equilibrium 
behavior even for slowly varying B and good thermal 
contact. These have been discussed extensively in con- 
nection with macroscopic quantum tunnelling. A well 
known example is Mni 2 -acetatei 2 i 22 i 23 i 24 i 25 

Departures from thermal equilibrium with the reser- 
voir due to imperfect thermal contact were observed 
and discussed for a wide range of sweep rates, dB/dt, 
from typical sweep rates in "dc magnets" to the very 
fast rates in pulsed magnets, 26 . 27 ! 28 . 29 ! 30 . 31 ! 32 ! 33 ! 34 The 
extreme case of a sample isolated from the thermal reser- 
voir (adiabatic conditions) was discussed by Wolf long 
ago, assuming thermal equilibrium within the sampled 
The magnetocaloric effect leads to cooling when any one 
of the MST's is approached. Such "cooling by adiabatic 
magnetization" has been observed many times in both 
dc fields and pulsed fields. 26 ! 29 . 32 ! 33 ! 36 



figure are the values of n for the finite chains responsible 
for each MST. 



F. Random Distribution 

Some of the preceding results used the probabilities 
P n for a random distribution of the magnetic ions. Al- 
though a random distribution is often found in diluted 
magnetic materials, non-random distributions were also 
observed. 1 . In the case of TMMC:Cd, difficulties of ob- 
taining a uniform Mn distribution were reported in the 
literatures* 1 ^ Therefore, the possibility of a non-random 
distribution cannot be ignored. The effects on the mag- 
netization curve caused by departures from a random 
distribution were discussed in Ref. For example, a 
tendency of the magnetic ions to cluster together is ex- 
pected to decrease M s . 



B. Some non-equilibrium models for pulsed 
magnetic fields 

1. Classification of non- equilibrium behaviors 

Non-equilibrium behavior in pulsed fields of millisec- 
onds duration is relevant to the interpretation of the 
present pulsed field data. For this purpose it will be use- 
ful to distinguish between three types of non-equilibrium 
situations: 

1. The spin-lattice relaxation is fast enough so that 
thermal equilibrium within the sample is estab- 
lished in a time which is very short compared to 
the pulse duration. The non-equilibrium behavior 
is then due to an inadequate heat flow between the 
liquid-helium bath and the sample, i.e., the sample- 
bath equilibration time is not short compared to 
the pulse duration. The "Inadequate Heat Flow" 
models are appropriate for this situation. 
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2. The time for establishing thermal equilibrium 
within the sample is not short compared to the 
pulse duration, but there is adequate sample-to- 
bath heat flow. The non-equilibrium is then gov- 
erned by the slow spin-lattice relaxation processes. 
There are many such processes. Here, only spin- 
flips associated with level crossings or with cross- 
relaxation (CR) will be discussedi 3 ^ 3 ^ 

3. Both the time for reaching equilibrium within the 
sample, and the time for reaching a sample-to-bath 
equilibrium are not short compared to the pulse 
duration. Some such cases were considered in the 
literature. 30 ! 3 ^ 33 ! 39 



2. Inadequate Heat Flow (IHF) Models 

The phonon bottleneck (PB) phenomenon has been 
known for a long time42*^i Treatments of this phe- 
nomenon in the context of the magnetization process of 
magnetic clusters led to the development of several IHF 
modelsi 27 i 29 i 30 i 31 i 33 i 34 The common feature of these mod- 
els is that the spin and phonon subsystems, within the 
sample, are very nearly in thermal equilibrium. The spin 
temperature T s , the phonon temperature, and the sample 
temperature are the same. However, the sample and the 
helium bath are not in thermal equilibrium (T s =/= Tbath)- 
Different IHF models treat the sample-to-bath heat flow 
differently. For a solid sample in contact with a liquid- 
helium bath the heat flow was assumed to be limited 
by the Kapitza resistancei 2 ^ 2 ^* 3 ^ Many of the qualitative 
features of the magnetization curve are common to all 
IHF models; they do not depend on the detailed treat- 
ment of the heat flow. 

Because the sample is internally in equilibrium, the 
equilibrium theory in Sec. [n] still applies. The absence 
of equilibrium with the helium bath enters only in the 
time dependence (and, hence, the -B-dependence) of the 
spin temperature T s . This dependence has been calcu- 
lated using various IHF models ^S^S^L^l Approach- 
ing the energy-level crossing (strictly, anticrossing) as- 
sociated with each MST, T s decreases. After passing 
through the anticrossing region, T s increases. If some 
heat flows between the sample and the bath then the lat- 
ter increase is large enough that T s is temporarily above 
Tbath- Miyashita and co- workers^ have called this be- 
havior of T s "the magnetic Foehn effect." 

The B-dependence T s can lead to the following qual- 
itative effects: 1) The dM/dB peaks are narrower than 
the thermal width. 2) The dM/dB peaks are asymmet- 
ric, i.e., the rise of dM/dB as the MST is approached is 
faster than the fall after passing through the MST. This 
is true both for increasing and decreasing B. 3) The 
magnetization and dM/dB exhibit hysteresis. 4) Under 
some conditions, a small "satellite" MST appears after 
the main MST. All these effects have been observed ex- 
perimentally. 



3. Cross Relaxation (CR) Model 

A severe non-equilibrium behavior, not explainable by 
IHF models, was observed in pulsed field experiments by 
Ajiro et alml It was interpreted in terms of CR between 
pairs and singles, and also between the pairs themselves. 
More recently, CR was discussed in the context of tunnel- 
ing theory? 38 ' 39 CR can involve both ground and excited 
states. 

CR is one of the mechanisms of spin relaxation. It 
involves simultaneous spin flips in weakly coupled clus- 
ters. The model of Ajiro et al. also includes a single spin 
flip in only one of the clusters. The later is actually not 
a CR process, and is better described as tunneling. In 
this model the spin relaxation rate is appreciable only at 
some values of B. 

The simplified picture used in Ref. |3J] ignored small 
level repulsions near level crossings. (Level repulsion 
is included in more detailed models that are based on 
tunneling. 38 39 ) In this simplified picture the relaxation 
rate can be fast only if simultaneous CR spin flips, or a 
single spin flip, do not change the total energy of the spin 
system. For a single spin flip, in one cluster, this hap- 
pens at energy-level crossings for this cluster. These level 
crossings include those of excited states. For CR between 
two weakly coupled clusters the energy-level separation in 
one cluster should match a level separation in the other. 
The two clusters may be of the same type (e.g., two cou- 
pled pairs) or of different types (e.g., a pair and a single). 
A similar criterion for energy-level separations applies to 
CR between three, or more, coupled clusters. 

In the equilibrium theory, the dM/dB peaks associated 
with MST's from Mn 2+ pairs occur at fields 

g/j, B Bi = 2i\J\, (15) 

where i=l, 2, . . . , (2S'-|-I). In the model of Ajiro et al. 
these peaks from pairs are called the "fundamental" 
peaks, and they are labeled as Pi. Specializing to pairs 
composed of Mn 2+ ions, there are five such peaks. In ad- 
dition to the fundamental peaks, other peaks are also pre- 
dicted. The most pronounced are the "second harmonic" 
peaks P m /2 at fields B m —m(Bi/2), with m=l to 10. 
The five peaks for odd m were observed clearly by Ajiro 
et al. The peaks for even m coincide with the fundamen- 
tal peaks. 

In addition to the second- harmonic peaks, fifteen third 
harmonic peaks Pfe/3 at Bk=k(Bi/3), were predicted. 
Many of these were also observed. Fourth and sixth har- 
monic peaks were also discussed. These results are from 
a model which involves only singles and pairs. The model 
can be extended to include other clusters. 

Some of the simultaneous spin-flip transitions which 
can give rise to P1/2 (the first of the ten second-harmonic 
peaks) are shown in Figs.QJa) and^Jb). Note that in ei- 
ther case, each of the two simultaneous spin flips increase 
the component of the spin along B by one unit. Figure 
shows a few level crossings that may contribute to var- 
ious peaks: the first three fundamental peaks Pi, P2, 
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FIG. 4: Some CR spin-flip transitions which may lead to 
the P1/2 peak in dM/dB. This peak is at B—\J\/gfiB- a) 
Simultaneous spin flips in two pairs, b) Simultaneous spin 
flips in a pair and in a single. 
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FIG. 5: Some energy-level crossings which may contribute 
to the "fundamental" peaks Pi, P2, and P3, and to the "har- 
monic" peaks P3/2 , P5/2, an d P2/3- The ¥ k / r peak is the 
fc-th peak of the r-th harmonic in the model of Ajiro et aim— 



P3; the third and fifth second-harmonic peaks, P3/2 and 
P 5 / 2 ; and the second third-harmonic peak, P2/3- Cross 
relaxation processes may also contribute to some of the 
same peaks, e.g., to P2/3- 



IV. EXPERIMENTAL TECHNIQUES 

A. Sample preparation 

The procedure of growing crystals of TMMC:Cd fol- 
lowed Ref. Il'ii The samples were grown by evapora- 
tion from water solutions of MnCl 2 .4H 2 0, (CH 3 ) 4 NC1, 
CdCl 2 .2H20, and some HC1 acid. The solutions were 
maintained at 30°C. As already noted, the Cd concentra- 
tion in the crystallized samples is about 50 times larger 
than in the starting solutioniiii^ Because the Cd con- 
centration in the solution decreases rapidly as the crys- 
tal growth progresses, a large solution volume (300 ml) 
was used to crystalize a "product" of TMMC:Cd with a 
total mass of about 100 mg. The large starting volume 
increased the uniformity of the Cd concentration in the 
product. 

Physically, the product consisted of many needle- 
shaped crystals. The long dimension of a needle (up to 
4 mm) was parallel to the hexagonal axis. The color grad- 
ually changed from pink towards white as the Mn con- 
centration x decreased. Each of the magnetization mea- 
surements, in both dc magnetic fields and in pulsed fields, 
used only a portion of the product, typically 30 mg. 

The Mn concentration x, was determined from the 
high-temperature dc susceptibility, x = M/B, measured 
using a SQUID magnetometer. The magnetic field B = 
0.1 T was well within the range where \ is independent 
of B. In the temperature range from 150 K to 300 K, the 
data for x were well described by a sum of a Curie- Weiss 
susceptibility and a constant due to the diamagnetism of 
the lattice. The concentration x was obtained from the 
Curie constant C. Strictly, the Curie- Weiss law is ac- 
curate only in the limit of very high temperatures. The 
percentage error in the Curie constant, resulting from the 
use of data between 150 K and 300 K, depends on the 
Mn concentration. Based on the results in Ref. l42l for 
the samples used here the error in C, and hence in x, was 
less 2%. The percentage error in the Curie- Weiss tem- 
perature 9 is larger than for C, but 9 does not enter into 
the determination of x. 

Values of x for several portions of the ~ 100 mg prod- 
uct obtained from a single solution were close to each 
other. A quantitative comparison between experimental 
dc magnetization data and theoretical simulations was 
carried only for samples 2, 4, and 5 in Tabled To increase 
the confidence in this comparison, values of x were deter- 
mined for the very same three samples. As a check, the 
Mn and Cd weight percents for these three samples were 
also determined directly by atomic emission spectroscopy 
with inductively coupled plasma (ICP-AES). The values 
of x deduced from ICP-AES were in reasonable agree- 
ment with the values from the susceptibility (see Tabled . 

The magnetization of several other samples (1*, 3*, 4*, 
and 5* in Tabled) was also measured in dc fields. How- 
ever, for these samples each value of x is from suscep- 
tibility data on a different portion of the same product. 
Samples 4 and 4* are two different portions of the same 
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TABLE I: Properties of the various samples. The Mn con- 
centration x, as determined from the magnetic susceptibil- 
ity (Suscept.) and from atomic emission spectroscopy with 
inductively coupled plasma (ICP-AES). The magnetic field 
Bi is at the first MST from pairs, and the magnetic field 
-Biquart is at the first MST from quartets. Both fields were 
determined from dc magnetization data. The NN exchange 
constant J was obtained from Bi. 



Sample 


X 


X 


Bi 


J/k B 


Biquart 


No. 


Suscept. 


ICP-AES 


(T) 


(K) 


(T) 


1* 


0.16 




8.85 


-5.94 




2 


0.22 


0.22 


8.70 


-5.85 


3.8 


2* 


0.22 










3* 


0.27 




8.85 


-5.94 




4 


0.25 


0.30 


9.10 


-6.11 


4.3 


4* 


0.28 




8.90 


-5.99 


4.5 


5 


0.48 


0.50 


9.65 


-6.48 


4.4 


5* 


0.50 




9.65 


-6.48 


4.4 



product, as are samples 5 and 5*. 

X-ray powder diffraction patterns were obtained at 
room temperature using Cu-K Q radiation. Data were 
taken on two samples from the same products as those 
of samples 1* and 4 (or 4*). The diffraction patterns for 
both samples were very similar to the pattern obtained, 
with the same equipment, on pure TMCC (x=0). No 
additional, or missing, diffraction peaks were observed. 
These results are consistent with a single crystallographic 
phase. 



B. Magnetization in dc magnetic fields 

Magnetization data in slowly varying magnetic fields 
(so-called "dc fields" ) were taken with a vibrating sample 
magnetometer (VSM). The VSM operated in 18 T super- 
conducting magnets. The sample was in direct contact 
with a liquid 3 He bath, which was in an insert dewar. 
The temperature 0.55 K was reached by pumping on the 
3 He bath. The field-sweep time (zero to 18 T) was about 
1 hour. 



C. Differential Susceptibility in pulsed magnetic 
fields 

The differential susceptibility, dM/dB, was measured 
in pulsed magnetic fields up to 50 T (500 kG). The 
techniques have been described earlieriSi The shape of 
the field pulse (B versus time) was approximately a half 
cycle of a weakly damped sine wave, with a rise time 
of 3.1 ms, and fall time of 4.3 ms. For each sample, 
data were first taken with the sample in the pickup coils, 
and shortly thereafter with the sample outside the pickup 
coils. The signal from the sample was obtained by taking 
the difference. 



The powder samples used in the pulse field experiments 
were obtained by crushing the mm-size needles of the 
growth products. Each sample, consisting of 20 to 30 mg 
fine powder (grain size less than 0.1 mm), was placed 
in a thin-walled (0.25 mm) cylindrical capsule made of 
Delrin. The capsule was immersed in a liquid 4 He bath 
which was maintained at 1.5 K. The Delrin capsule had 
a small hole at its bottom. The hole (covered by a tissue 
paper) allowed a direct contact between the sample and 
the bath of supcrfluid helium. However, previous exper- 
iments have indicated that despite such a direct contact 
the sample may not be in thermal equilibrium with the 
bath during the 7.4 ms field pulseiSS* 3 ^ 



V. MAGNETIZATION IN DC MAGNETIC 
FIELDS 

A. Experimental Results 

1. Gross Features 

Figure Ela) shows magnetization data at 0.55 K for 
samples 2, 4, and 5. These "dc data" were actually ob- 
tained with a sweep rate of ~ 0.3 T/min. No hysteresis 
was observed. The (very small) corrections for lattice dia- 
magnetism and addenda are included in Fig. so that 
the magnetization M is that of the Mn 2+ ions. The nu- 
merical derivatives, dM/dB, of these curves are shown 
in Fig. IHIb). The Mn concentrations x, measured on the 
very same samples, are given in Tabic [I] 

The main features in Fig. E^a) are: 1) a fast rise of 
M at low H, 2) a large MST near 9 or 10 T, and 3) 
a smaller MST near 4 T. There are also indications of 
other MST's at higher fields. For example, sample 2, 
for which the data extend to slightly higher fields than 
for the other samples, shows the beginning of a large 
MST near the top of the field. These main features of 
the dc magnetization curves agree with theoretical pre- 
dictions for the equilibrium magnetization, such as those 
in Fig. El^b) . The range of the reduced field b that cor- 
responds to the experimental data in Fig. [5] extends up 
to about 4. The exact maximum value of b is slightly 
different for different curves in this figure. 

The fast rise of the M at low B corresponds to the 
alignment of the zero-field-ground-state spin, S*t(0), of 
finite chains with odd n. The main contribution to this 
fast rise is from the singles (n=l). The large MST ob- 
served near 9 or 10 T corresponds, essentially, to the first 
MST from pairs (clusters with n=2). Other contributions 
to this observed MST are from some longer chains that 
have a MST at nearly the same field. For example, the 
second MST from quartets (chains with n=4) is predicted 
to occur at a field which is only 2% higher than that of 
the first MST from pairs. Neither the triplets {n=3) nor 
the quintets (n=5) have a MST near this field. The total 
contribution of chains with n>2 to the observed MST 
near 9 or 10 T is expected to be smaller than the contri- 
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third MST from the quartets is also predicted to occur 
within the field range of Fig. El However, the predictions 
for this MST were not fully confirmed by the data. For 
example, in Fig. El[b) the derivative curve for sample 5 
exhibits a small peak near 15 T. For the same sample, 
the predicted field at the third MST from the quartets is 
higher by about 1 T. 

Some features of the experimental results in Fig. de- 
pend on the Mn concentration x. The first MST from 
the quartets stands out more clearly as x increases. This 
trend is expected from the probability curves in Fig. ^ 
The cluster populations N n are related to these proba- 
bilities by Eq. @. As x increases there is an increase 
in population ratios N4/N2 between quartets and pairs, 
and N4/N1 between quartets and singles. Therefore, as x 
increases, the MST from quartets stands out more clearly 
in comparison with the MST from pairs, and also in com- 
parison with fast magnetization rise at low fields. 

Another feature that depends on x is the value of the 
magnetic field at the large MST near 9 or 10 T. In 
Fig. HJb) the peak associated with this MST shifts to 
slightly higher fields as x increases. The field at this peak 
is expected to be very close to B\. Numerical values as 
a function of x are listed in Table D The change of B 1 
is attributed to a slight dependence of the NN exchange 
constant J on the Mn concentration. 



FIG. 6: (a) Magnetization M of samples 2 (x=0.22), 4 
(x=0.25), and 5 (x=0.48), measured at 0.55 K in dc magnetic 
fields. The results have been corrected for lattice diamag- 
netism and addenda. The SI unit Am 2 /kg is equivalent to 
1 emu/g. (b) The numerical derivative, dM/dB of the mag- 
netization traces. The derivative curves for samples 4 and 5 
have been shifted upwards. 



bution from the pairs. The reason is that for x=0.5 the 
populations N n of these longer chains are small compared 
to the population N2 of the pairs. 

For sample 2, a substantial portion of the second MST 
from pairs is also seen at the highest fields. The peak 
in the derivative dM/dB for this sample, near the top of 
the field in Fig.EJb), is close to the expected field for this 
MST, i.e., 52=251 = 17.4 T. The beginning of the second 
MST from pairs is also seen in the derivative curve for 
sample 4. 

The relatively small MST near 4 T, seen in Figs. Efa) 
and (b), is identified as the first MST from the chains 
with n=4 (so-called "string quartets"). The first MST 
from the quartets is predicted to occur at B—0A75Bi, 
where B\ is the field at the first MST from pairs. The 
experimental results in Table [i] are in reasonable agree- 
ment with this prediction. The second MST from quar- 
tets which, as already noted, is predicted to occur at a 
field which is only 2% higher than B\, was not resolved 
at 0.55 K. This was expected because at 0.55 K the 
broadening of any MST due to the finite temperature 
("thermal broadening") is more than 10% of B\. The 



2. NN Exchange Constant 

The NN intrachain exchange constant J was obtained 
from the field B\ of the first MST from pairs, using 
Eq. jTSjl and assuming g=2.00 for the Mn 2+ ion. Val- 
ues of J for all the samples are given in Table [I] The 
exchange constant for x=0.5, J/k^ — — 6.5 K, is about 
10% higher than for x=0.22. This 10% change is too 
large to be accounted for by the unresolved MST from 
the quartets. As already mentioned, J/k^ = —6.6 K for 
pure TMMC (x=l). We are not aware of any theoretical 
calculation of the x-dependence of J in this system. 



B. Comparison with Simulations 

Figures [7HH1 compare the experimental results with 
numerical simulations based on equilibrium theory for 
the single- J model ( Section ITXp . The comparison is for 
the normalized magnetization m=M/ Mq, where Mq is 
the true saturation magnetization. The "experimental" 
curves use the measured M (Fig. and the calculated 
saturation magnetization Mq for the Mn concentration 
x. The simulations assume a random distribution of the 
Mn ions, and use the values of J and x(Suscept.) given in 
Table |U for that sample. Clusters with n < 5 are treated 
exactly, and the rise-and-ramp approximation is used for 
the total contribution of larger clusters. There are no 
adjustable parameters in the simulations. 
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FIG. 7: Comparison between the dc data for x=0.22 at 0.55 K 
with a simulation based on equilibrium theory for the single- J 
model. The left ordinate scale is for the normalized magneti- 
zation m=M /Mo, where Mo is the true saturation value. The 
right ordinate scale is for dm/dB. Solid curves are from the 
data in Fig. [5] and the calculated Mo- The dotted curves are 
from the simulation. Only thermal broadening, at the actual 
temperature 0.55 K, is included in the simulation. 



For sample 2, with x=0.22, the agreement between ex- 
periment and theory is reasonably good (Fig. 0). The 
difference between the measured and simulated magne- 
tizations is a few percent. It is comparable to the total 
experimental uncertainty, mainly from the uncertainty 
in x. The observed magnitude of the first MST from 
the pairs is in agreement with the simulation. How- 
ever, the associated experimental dM/dB peak is some- 
what broader than in the simulation. Physical mecha- 
nisms that broaden MST's were discussed in Ref. [lj. In 
the present case thermal broadening at the experimental 
temperature, T=0.55 K, is expected to be the strongest 
of these mechanisms. It is the only broadening mecha- 
nism that was included in the simulations. Non-thermal 
causes of line broadening include the dipole-dipole in- 
teraction and local strains associated with the random 
replacement of Mn by Cd. Because non-thermal broad- 
ening was neglected in the simulation, it is not surpris- 
ing that the experimental dM/dB peaks are somewhat 
broader. The numerical differentiation of M with respect 
to B also broadens the experimental peak slightly. 

Figure |S| shows that for x=0.25 the agreement between 
experiment and simulation is, again, reasonably good. 
However, the observed dM/dB peak at 4.3 T, from the 
first MST for quartets, is somewhat larger than expected. 
The simplest interpretation is that the number of quar- 
tets is larger than given by a random distribution. The 
behavior of the derivative dm/dB near 15 T is attributed 
to the third MST from quartets. 

Figure[5]shows that the agreement between experiment 



FIG. 8: Comparison between the dc data for x=0.25 at 0.55 K 
(solid curves) with a simulation based on equilibrium theory 
for the single-J model (dotted curves). The ordinate scales 
are for m—M/Mo, and dm/dB. 
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FIG. 9: Comparison between the dc data for x=0.48 at 0.55 K 
(solid curves) with a simulation based on equilibrium theory 
for the single-J model (dotted curves). The ordinate scales 
are for m=M/Mo, and dm/dB. 



and simulation for x=0.48 is only fair. The measured 
magnetization is somewhat smaller than expected, over 
the entire field range. In particular, the initial rise of the 
measured magnetization is smaller than in the simula- 
tion. This discrepancy suggests that the Mn ions have a 
slight tendency to bunch together, compared to a random 
distribution. 1 Also, the peak of dm/dB near 15 T is at 
a lower field than predicted for the third MST from the 
quartets. The reason for this behavior is unclear. In our 
view, despite these discrepancies the agreement between 
experiment and theory is still reasonable. 
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VI. DIFFERENTIAL SUSCEPTIBILITY IN 
PULSED FIELDS 

Pulsed- field data for dM/dB were taken on powder 
samples with x w 0.50 (from the same product as sam- 
ples 5 and 5*), x w 0.22 (same product as 2 and 2*), and 
x « 0.16 (same product as 1*). As noted in Sec. HVI each 
sample was in direct contact with a liquid-helium bath, 
maintained at Tb a th=l-5 K. However, such a direct con- 
tact does not ensure thermal equilibrium with the bath 
during the 7.4 ms pulse. 

The pulsed- field data are presented next. Non- 
equilibrium features are pointed out, and many of them 
are then interpreted. However, some aspects of the non- 
equilibrium behavior are still not fully understood. 

A. Experimental Results in Pulsed Fields 

Figure ITU1 shows some dM/dB data for x « 0.50. Part 
(a) of this figure shows an "up trace" (increasing B) and 
a "down trace" (decreasing B) obtained during the same 
field pulse. Part (b) gives an expanded view of a portion 
of the down trace. 

Results for x « 0.22 are shown in Fig. Illf a). This 
figure covers the wide range of dM/dB values that is 
required to display the large hysteresis in low fields. Data 
taken during the same field pulse, but which cover a much 
narrower range of dM / dB values, are shown in Fig JllT b). 

Results obtained for x « 0.16 during a pulse with 
a maximum magnetic field B max ~l7 T are shown in 
Fig. I12f a) . Figure I12f b) shows the decreasing-field por- 
tion of a trace obtained during another field pulse for 
which the maximum field was 34 T. [Unlike all the 
other pulsed-field data in the present paper, the data 
in Fig. I12r b) have not been corrected for the monotonic 
background, because no "background shot" (with the 
sample out of the pickup coils) was taken in this case.] 
Figure 1131 shows an expanded, and slightly smoother, 
view of the field-down portion of the Fig. I12f a). 

B. Discussion of non-equilibrium effects 

The pulsed-field data do not provide any new informa- 
tion about the exchange constant J. However, these data 
show interesting non-equilibrium effects. The discussion 
below focuses primarily on: 1) those features of the data 
that indicate the absence of thermal equilibrium; 2) the 
change of the non-equilibrium behavior with the Mn con- 
centration x; and 3) physical mechanisms that can give 
rise to such non-equilibrium effects. 

1. x«0.50 

Consider first those pulsed-field results in Fig. II Ul that 
are below 15 T. There is a prominent peak near 10 T. It 



stands out more clearly in the down trace (decreasing B) 
than in the up trace. This peak corresponds to the first 
MST from pairs. The down trace also shows a smaller 
peak just below 5 T. This smaller peak corresponds to 
the first MST from the quartets. Although both of these 
peaks were also observed in the dc data for the equi- 
librium magnetization [Fig. E[b)], at least two features 
of the pulsed-field data indicate departures from equilib- 
rium with the helium bath, at Tb a th=l-5 K. First, the up 
and down traces below 15 T are different, i.e., there is an 
hysteresis in this field range. Second, the widths of both 
the peak near 10 T and near 5 T are substantially smaller 
than the equilibrium width at the bath temperature. 

The "width" of a MST will always refer to the full 
width at half maximum of the associated peak in dM/ dB. 
The various contributions to the equilibrium width were 
discussed in Ref. QJ. Often, temperature broadening is 
important. The thermal width at the temperature T is 

{5B) T = 3.53fc B T/. gMB . (16) 

Because non-thermal broadening mechanisms are also 
present in equilibrium, the thermal width is a lower limit 
for the actual equilibrium width. The calculated ther- 
mal width, (SB)t=3.9 T for the actual bath temperature 
7bath=l-5 K, is shown in Fig.HuTb). 

To obtain the experimental width it is necessary to 
choose a baseline for the peak in dM/dB. Such a choice 
is not always obvious. In Fig. llOf b). two possible choices 
for the peak near 10 T are indicated by the dashed lines 1 
and 2. Baseline 1 leads to a full width at half maximum of 
2.2 T. Baseline 2 leads to a width of 1.9 T. On this basis 
we conclude that the experimental width is substantially 
smaller than the equilibrium width. 

The choice of a baseline for the peak near 5 T is also 
not obvious. The particular choice shown as a dashed 
line leads to a width of 1.2 T. Although other choices 
may lead to a larger experimental width, it seems that 
any reasonable choice will lead to a width that is smaller 
than the thermal width. Thus, the experimental widths 
of both of the peaks in Fig. HUf b) are smaller than the 
width that would have occurred had the sample been in 
equilibrium with the helium bath. 

Figure HUI also shows the second MST from pairs, near 
20 T. This peak is broader than the first peak from 
the pairs. A significant feature of this peak is that it is 
asymmetric. As a function of time, the rise of dM/dB as 
the peak is approached is faster than the fall after passing 
through the peak. This asymmetry is observed in both 
increasing and decreasing B. An asymmetry of this type 
is expected from models that assume an inadequate heat 
flow (Sec. IIIIBII . The data in Fig. ^3 therefore suggest 
that the non-equilibrium behavior for x w 0.50 is due to 
inadequate sample-to-bath heat flow. 

The IHF scenario can also account for the hysteresis 
below about 15 T. When the thermal contact with the 
bath is poor, strong magneto-caloric effects are expected 
from the singles (clusters with n=l). These effects are 
analogous to those involved in the cooling of a paramag- 
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net by adiabatic demagnetization (and warming by adia- 
batic magnetization). Of course, the actual processes in 
the present case are not truly adiabatic, because there is 
some heat flow between the sample and the bath. In our 
view, no feature of the data in Fig.llOlreauires CR for its 
explanation. The observed non-equilibrium effects seem 
to be explainable by an inadequate sample-to-bath heat 
flow. 

Another issue (not directly related to the non- 
equilibrium behavior) involves the third MST from pairs, 
seen in Fig. linear 30 T. Compared to the first and sec- 
ond MST's from pairs, the third MST is less well defined. 
Specifically, dM/dB hardly decreases on the high- field 
side of the third MST. This behavior is explained by the 
three small MST's from quartets, triplets, and quintets 
that are expected between the third and fourth MST's 
from pairs (see the MST's near 6=7 in Fig.[3J). For x=0.5 
the predicted combined size of these three small MST's 
is comparable to the size of one MST from pairs. 



2. N on- equilibrium behavior for x ~ 0.16 

The absence of thermal equilibrium with the helium 
bath is very evident in the pulsed-field data for x « 0.16. 

(i) A pronounced hysteresis is seen in Fig. I12f ah 

(ii) In the down portion of the pulse (see Fig. 1 1 3|> the 
widths of the peaks near 10, 5, and 2.7 T are 1.4, 0.9, 
and 0.6 T, respectively. These values are small com- 
pared to a thermal width of 3.9 T at Tb a th=l-5 K. 
A width that is smaller than the thermal width im- 
plies a non-equilibrium behavior. 

(iii) The peak near 5 T is very pronounced in the pulsed- 
field data shown in Figs. ll2f a1 and 1 131 but is barely 
seen in dc data on a sample from the same product. 
The dc data, shown in Fig. should be represen- 
tative of equilibrium behavior. Therefore, the pro- 
nounced peak in the pulsed field data is regarded as 
a non-equilibrium effect. 

(iv) The second MST from pairs, near 18 T, is barely ob- 
served in Fig.UUb), and the third MST from pairs, 
near 27 T, is totally absent. The corresponding 
field-up trace (not shown) exhibits a similar behav- 
ior. Once again the pulsed field data are contrasted 
with the equilibrium magnetization data in Fig. 1141 
The latter data, which extend up to 17.5 T, show 
a significant portion of the second MST from pairs, 
and they indicate that in equilibrium the sizes of 
the second and first MST's from pairs are compa- 
rable. That is, unlike the behavior in pulsed fields, 
the second MST is not small compared to the first. 
In equilibrium, all MST's from pairs are expected 
to be comparable, which is inconsistent with the re- 
sults of Fig.H2Tbj. 
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FIG. 10: (a) Pulsed-field data of dM/dB for x « 0.50 in both 
increasing B ("up") and decreasing B ("down"). The up and 
down traces are shifted vertically relative to each other, (b) 
Expanded view of a portion of the down trace. The calculated 
thermal width (SB)t at the bath temperature Xbath=l-5 K is 
indicated. Dashed lines show some choices of baselines used 
to obtain the experimental widths of two peaks. 



3. IHF and CR Scenarios for x « 0.16 

The large hysteresis at low fields [Fig. I12f a1] is not 
uncommon in pulsed field experiments on diluted mag- 
netic materials M^ 7 - In the present case the hysteresis is 
largely due to the slow response of the singles to the 
rapidly changing B. In the up portion of the pulse the 
alignment of the spins of the singles is not completed un- 
til B is above 10 T. In the down portion of the pulse 
these spins remain largely aligned until B drops below 
1 T. For x w 0.16, the majority of the spins (70% for a 
random distribution) are singles, so that the hysteresis is 
very pronounced. 

MST's from pairs and larger clusters stand out more 
clearly in the down portion of the field pulse because the 
singles remain largely aligned until B drops below 1 T. 
For this reason the field-down portion of the traces in 
Figs. 1131 and I12f b) is discussed first. 
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FIG. 11: Pulsed-field data of dM/dB for x « 0.22. (a) Over- 
all view, showing the large hysteresis below about 10 T. (b) 
Expanded view of the field-up and field-down traces obtained 
during the same field-pulse. The low-field portion of the up 
trace is excluded. The up and down traces are not shifted 
vertically relative to each other. 



FIG. 12: (a) Results for x ~ 0.16, obtained during a field 
pulse with a maximum field -B m ax = 17 T. The up and down 
traces are not shifted vertically relative to each other, (b) The 
down portion of a trace for x ~ 0. 16, obtained in another pulse 
with -B max =34 T. This particular trace, unlike all others, is 
not corrected for background. 



Spin-lattice relaxation times are often longer for lower 
x, which is the likely cause of the more pronounced non- 
equilibrium behavior for x « 0.16 compared to x w 0.50. 
An important issue in the data interpretation is whether 
the spin-lattice relaxation is fast enough to maintain 
equilibrium within the sample. In that case the IHF sce- 
nario would apply (see Sec. IIII B 2jl . An alternative is 
a more severe non-equilibrium behavior which is better 
described by the CR scenario, including single spin flips 
near level crossing s of excited states ('Sec. IIII B 

Below we consider both the IHF and CR scenarios 
for the down portion of the pulse. The preferred sce- 
nario cannot be chosen on the basis of the fields at which 
the MST's occur, because the differences are often small 
compared to the experimental accuracy. However, the 
two scenarios lead to different relative sizes of the peaks 
in dM/dB. The poor agreement of the observed rela- 
tive sizes with those predicted by the IHF scenario will 
suggest that the CR scenario is preferable for this Mn 
concentration. 



Consider first the IHF scenario. In this scenario the 
magnetization of each cluster type is the equilibrium 
magnetization at T s . For the low T s indicated by the 
small widths of the observed MST's (in the down portion 
of the pulse), this magnetization is that of the ground 
state. The largest peak in Fig. ^1 just below 10 T, is 
mainly due to the first MST from pairs. The second 
largest peak, near 5 T is the first MST from quartets. 
It is predicted to occur at a field which is 0.48 times 
that of the first peak for pairs. The second dM/dB peak 
from quartets should be at a field which is 2% higher 
than the first peak from pairs. The structure of the peak 
near 10 T is possibly due to a superposition of these two 
MST's, although it is much wider than 2%. 

In the IHF scenario the field at the first MST for finite 
chains with even n can be estimated from the energy 
levels given in Ref. [23- For n between 4 and 10, this field 
is nearly proportional to l/n. A similar approximate 
dependence on n holds for AF rings (closed chains) £ On 
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FIG. 13: An expanded, and slightly smoother, view of the 
field-down portion of Fig. I12f a~l . 
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FIG. 14: The dc magnetization M at 0.55 K for x w 0.16 
(sample 1*). Also shown is the numerical derivative dm/dB. 
In the text, some features of these equilibrium-magnetization 
data are contrasted with pulsed-field data on a similar sample 
(Figs. [H and ESJ. 



this basis the small sharp peak near 2.7 T (in Fig. I12|l is 
due to octets or sextets. There is also an indication of 
a small peak near 2.0 T, which would be attributed to 
larger clusters. 

Although the IHF scenario accounts for the fields of 
many of the observed MST's, this scenario is very ques- 
tionable for this Mn concentration. Assuming that the 
Mn cations are randomly distributed, the number of 
quartets for x = 0.16 is smaller than the number of pairs 
by a factor of 39. For sextets and octets the factors are 
1.5 x 10 3 and 6 x 10 4 , respectively. Therefore, unless 
the deviations from random distribution are extremely 



large, it should not have been possible to observe MST's 
from sextets or octets if the behavior followed the IHF 
scenario. The first MST from quartets might have been 
detectable, but it should have been very small compared 
to the MST from pairs. This was not the case in the 
pulse field experiments (Figs. IT21 and 113( 1 . As predicted, 
the equilibrium data in Fig. 1141 for nearly the same x, 
indicate that the first MST from quartets is much less 
pronounced than the first MST from pairs. 

In Fig. 112( b) the first MST from pairs stands out 
clearly but the second MST from pairs is barely visible, 
and the third is totally absent. These results also are 
not well understood within the IHF scenario, although 
some indication of such a behavior appeared in simula- 
tions by Nakano and Miyashita for iron clusters with ring 
structured 

In the CR scenario the MST's in Figs. H^b) and 1131 

(both for decreasing B) are interpreted as follows. The 
peak just below 10 T is the first fundamental peak Pi, 
with some contributions from P2/2, P3/3, etc. The peak 
at 5 T is the second harmonic peak Pi/2- The small 
peak at 2.7 T is the third-harmonic peak P1/3 or the 
fourth-harmonic peak P1/4. The broad peak near 6.6 T 
is possibly P2/3- Because this interpretation uses only 
pairs and singles*^ it is not open to objections based on 
the low populations of quartets and larger clusters. The 
CR scenario also accounts for some features of the up 
trace in Fig. 112( a). The large peak near 5 T is attributed 
to the CR process of the type shown in Fig.0|b), except 
that the directions of both spin flips are reversed. This 
process allows the singles to relax toward a state with 
a higher magnetization. The magnetization of the pairs 
also increases by this process. 

The process in Fig. 0Jb) involves only one single and 
one pair, and is the simplest CR process between singles 
and pairs. It accounts for the Pi/2- More complicated 
CR processes can lead to other "harmonic peaks," such as 
Pi/3i P2/37 P3/4- These peaks are expected to be smaller 
than Pi/2- The peak P1/3 may involve processes such as 
a CR between a pair and two singles, or between a pair 
and a single which undergoes a double spin flip. The peak 
P2/3 may involve a spin flip in a single and spin flips in 
two pairs, etc. In the up trace the observed peak at 3.2 T, 
and the small peaks at 6.4, and 7.4 T may correspond to 
the P1/3, P2/3 and P3/4 harmonics. The peak observed 
near 10 T is a superposition of the first fundamental peak 
from pairs, Pi, and the harmonics P2/21 P3/37 etc. 

As already noted, in Fig. I12f the second MST from 
pairs is barely visible, and the third is totally absent. 
These results suggest that for x « 0.16, the spin relax- 
ation for pairs in fields above 15 T is very slow compared 
to a millisecond. As B sweeps through a region where 
a MST from pairs should have occurred, the pairs are 
unable to relax towards the new ground state. We spec- 
ulate that the slow spin relaxation for pairs is mainly due 
to a reduction of CR between pairs and singles, and that 
this reduction is related to the saturation of the singles 
in fields above 15 T. CR between different pairs, or be- 
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tween pairs and larger clusters, is expected to become 
slower as x decreases^ Among the three samples, such 
CR processes should be least efficient for x » 0.16. 

4- N on- equilibrium behavior for x ~ 0.22 

Non-equilibrium behavior is also evident for x « 0.22. 
The large low-field hysteresis for this Mn concentration, 
in Fig. Illf a). is somewhat similar to the hysteresis in 
Fig. I12f a') for x « 0.16. In the up trace, the large peak 
near 4.5 T is identified as the P1/2 peak, and is attributed 
to the cross relaxation process in Fig.^Jb), with the ar- 
rows reversed. 

The down trace in Fig. Illf b) exhibits large peaks near 
8.9 T and 4.4 T approximately, and a small peak near 
2.1 T. The widths at half height of these peaks, 1.3, 0.9, 
and 0.7 T, respectively, are all much smaller than thermal 
width of 3.9 T at Tbath- These widths, which are similar 
to those for x « 0.16, indicate non-equilibrium behavior. 

The largest peak in Fig. lllf b). at 8.9 T, is undoubtedly 
the first fundamental peak from pairs, Pi. The second 
fundamental peak P2 is also observed near 18 T. The 
CR scenario predicts large second-harmonic peaks P1/2, 
P3/2, P5/2 at 4.5, 13.4, and 22.3 T, respectively. These 
fields are close to 4.4, 13.4, and 21.8 T, of large peaks in 
Fie. II If hi. The peak observed near 2.1 T is consistent 
with P1/4. It is possible, but far from certain, that the 
small peaks at 6.2 and 12.0 T are P2/3 and P4/3. 

The main difference between x w 0.22 and x « 0.16 
is that for the higher Mn concentration there are still 
prominent MST's above 10 T. For x « 0.50 the second 



and third MST's from pairs are even more pronounced 
[see Fig. ITUT a)]. These results suggest that at these fields 
the spin relaxation rate for pairs increases rapidly with 
Mn concentration. We tentatively attribute this trend 
to the expected increase with x of the efficacy of CR 
processes involving pairs, and pairs and larger clusters. 



5. Summary of the analysis of non- equilibrium behavior 

A definitive interpretation of the observed non- 
equilibrium behavior in pulsed fields is still lacking. How- 
ever, it appears that for x s=s 0.50 the non-equilibrium 
behavior is better explained by the IHF scenario. For 
x w 0.16 and x « 0.22 the data are better explained 
by spin flips associated with CR and with level cross- 
ings. The data suggest that at the high magnetic fields 
where the magnetization of singles is saturated, the spin 
relaxation rate for pairs increases rapidly with increas- 
ing x. This increase is tentatively attributed to the x- 
dependence of CR processes involving pairs. 
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